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We show theoretically that arbitrary coherent rotations can be performed quickly (with a gating 
time ~ 1 ns) and with high fidelity on the spin of a single confined electron using control of exchange 
only, without the need for spin-orbit coupling or ac fields. We expect that implementations of this 



scheme would achieve gate error rates on the order of 77 < 10 
of several known error-correction protocols. 



The elementary building-blocks for universal quantum 
computing are a two-qubit entangling operation, such 
as the CNOT-gate or \/swAP-gate and arbitrary single- 
qubit rotations. For qubits based on single electron spins 
confined to quantum dotsji recent experiments have now 
achieved the two-qubit VswAP-gate^ and single-spin co- 
herent rotations.— If these operations are to be used in a 
viable quantum information processor, they must be per- 
formed with a sufficiently small gate error per operation 
77 <C 1. The threshold values of 77 required for effec- 
tive quantum error correction depend somewhat on error 
models and the particular error-correction protocol, but 
current estimates are in the range 77 < 10~ 2 — 10~ 4 4^ To 
achieve these low error rates, new schemes must be devel- 
oped to perform quantum gates quickly and accurately 
within the relevant coherence times. 

Previous proposals^ and recent implementations^ for 
single-spin rotation have relied on ac magnetic fields to 
perform electron-spin resonance (ESR) . In ESR, difficul- 
ties with high-power ac fields limit single-spin Rabi fre- 
quencies to values that are much smaller than the op- 
eration rates typically associated with two-qubit gates 
mediated by exchanged To circumvent these problems 
while still achieving fast coherent single-qubit rota- 
tions, there have been several proposals to use ex- 
change or electric-field (rather than magnetic-field) con- 
trol of electron spin states. These proposals aim to 
perform rotations on multiple-spin encoded qubitSf^ 
or require strong spin-orbit interaction; 9 - ' 10 ' 11 ' 12 cou- 
pling to excited orbital states^ or rapid pulsing of 
magnetic fields. 14 Qubits encoded in two states hav- 
ing different orbital wave functions are susceptible to 
dephasing through fluctuations in the electric environ- 
ment, even in the idle state i 15 ' 16 ' 17 Proposals that make 
use of the spin-orbit interactio n 9 ' 10 ' 11 ' 12 are restricted 
to systems where the spin-orbit coupling is sufficiently 
strong, excluding promising architectures such as quan- 
tum dots made from Si:SiGei£ and carbon nanotubes 
or graphene sheets! 19 ' 20 ' 21 Sufficiently rapid pulsing of 
magnetic fields 14 may not be feasible in GaAs, where 
the electron-spin coherence time is on the order of t c ~ 
10nsA2£ 

Here we propose to perform single-qubit rotations in 
a way that would marry the benefits of demonstrated 
fast electrical control of the exchange interaction^ with 
the benefits of naturally long-lived single-electron spin 
qubitsi Our proposal would operate in the absence of 



in GaAs quantum dots, within reach 
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Figure 1: (Color online) Possible setup to implement the 
scheme proposed here. Ancillary electron spins at z\ are main- 
tained in a polarized state with a large Zeeman field b\ along 
z. Qubit spins at Z2 are free to precess in a weaker effective 
Zeeman field lying in the x-z plane: A = (62, 0, 6| — J/2). 
Here, J is the exchange coupling between the qubit and an- 
cillary spins and b2 is the qubit Zeeman field in the absence 
of exchange. When b\ 3> 6| 3> °2, z-rotations are performed 
if J ~ and ^-rotations are achieved when J ~ 26|. 



spin-orbit coupling and would act on single electron spins 
without the use of ac electromagnetic fields, in the pres- 
ence of a fixed Zeeman field configuration (Fig. [J). This 
scheme applies to confined electrons in any structure with 
a locally controllable potential. Specifically, this scheme 
may be applied to electrons above liquid helium, bound 
to gated phosphorus donors in silicon, and in quantum 
dots formed in a GaAs two-dimensional electron gas, 
nanowires, carbon nanotubes, or graphene. 

Hamiltonian-We begin from a standard tunneling 
model for the two lowest orbital levels of a double quan- 
tum dot, including tunnel coupling t±2, on-site repulsion 
U c , nearest- neighbor repulsion U' c , local electrostatic po- 
tentials Vi o\ and a local Zeeman field t>i(2) on dot 1(2) 
(see Refs. J23j2i[ and references therein): 



n = - Vim* + U *Y1 



In 



h2 J2 ( d ld2* + 4 a d la ) - £ S, . b, . (1) 



Here we have set H = 1, di a annihilates an electron in dot 
I = 1,2 with spin a, 7i/ CT = d, di a is the usual number 
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operator, and Si = | J2 S s' c i a tr s,a'Cls> is the spin density 
on dot /. We choose |e±56 z | > \t 12 \, \6b z \ > |i 12 |, with 
<56 z = (6f - 6|)/2 and e = Va - Vi - U c + U' c , which 
favors the (1,1) charge state (where {N\,No/) denotes a 
state with N^) electrons on dot 1(2), see Fig. |2J|. Ad- 
ditionally, we require a large Zeeman field along z in dot 
1 (\b\ | 3> |&i ,3/ |) so that the spin on dot 1 is frozen into 
its spin-up ground state. For simplicity, we furthermore 
choose b\ = 0. Eq. |T]) then reduces to the following 
low-energy effective Hamiltonian for the spin on dot 2: 



J(e)/2). (2) 



When |e| > \8b z \, J(e) « -2t\ 2 /e. Thus, for a fixed 
Zeeman field b 2 , the direction and magnitude of the ef- 
fective field A can be tuned with gate voltages via its 
dependence on e (see Fig. HJc)). Eq. ([2]) follows di- 
rectly from a much more general Hamiltonian of the 
form H = J(e)Si • S2 — Y]r hi ■ Si in the limit where 
|t>i| 3> |b2|,J, and so this scheme is not limited to 
the particular Hamiltonian given in Eq. ([1]), which ne- 
glects the long-ranged nature of the Coulomb interaction 
and excited orbital states. The long-ranged part of the 
Coulomb interaction (the exchange integral) contributes 
a small fraction to J(e) compared to the tunneling con- 
tribution when the out-of-plane magnetic field is zerc 25 
and contributions to J(e) due to excited orbital states^ 
are a small correction when |e| < J(o,2)-> where J(o,2) is 
the single-dot exchange coupling on dot 2. Outside of 
this range of validity, the functional form J(e) could be 
obtained empirically, as has been done in Ref. (26j . 

Qubit gates-Arbitrary single-qubit rotations can be 
achieved with the appropriate composition of the 
Hadamard gate (H) and it/ 8 gate (T)21: 



H = 



1 



1 

-1 



T = 







(3) 



Up to a global phase, T corresponds to a rotation about 
z by an angle 4> = tt/4. This operation can be performed 
with high fidelity by allowing the qubit spin to precess 
coherently for a switching time t s = </>/A z at the operat- 
ing point tA in Fig. Ufa), where A z ^> A x . The H gate 
can be implemented by pulsing e (see Fig. djc)) from ca, 



where A 2 w bo to ts 



-n 2 /b s 2, 



where A z and back. 



The pulse is achieved with a characteristic rise/fall time 
t, and returns to e = 6a after spending the pulse time 
t p at e — eb- If 6f <C TloS induces approximate z- 
rotations during the rise/fall time, and x- rotations when 
e = €b- The entire switching process (with total switch- 
ing time t s = t p + 4r) is described by a time evolution 
operator U = Te^o* 3 dtA(t)-a/2^ wmdlj for b x < ^ is 
thus approximately given by 



U « U(ct> x , 4> z ) = R z R x (-cf> x ) R z \-^.\ , (4) 

where <p x = A x t p and 4> z — J * s dtA z (t). Here, Rh(4>) is 
a rotation about the n-axis by angle (p. When <f> x = it/ 2 
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Figure 2: (Color online) (a) Charge stability diagram indicat- 
ing the ground-state charge configuration (Ni,N2) for local 
dot potentials V\, V%. In the (1, 1) configuration, the exchange 
interaction J(e) can be tuned by shifting the double-dot po- 
tential difference e ~ V% — V\. (b) When the electron spin in 
dot 1 is polarized, the qubit electron acquires a Zeeman shift 
given by tfa/e = — J(e)/2 due to virtual hopping processes 
that are allowed for spin-down, but forbidden for spin-up due 
to the Pauli principle, (c) Energy spectrum of the Hamilto- 
nian given in Eq. fTJ in the presence of a strong inhomoge- 
neous magnetic field. 



and <f> z = 7T, Eq. (g]) gives an H gate, up to a global 
phase: U(tt/2,it) = iH. 

Errors— We quantify gate errors with the error rate r\ = 
I — J 7 , where T is the average gate fidelity, defined by 



T = — I dftTr ( Up m (6, ct>)WUp in (0, 4>)W). (5) 



Here, p- u 



,(/>) {0,<f>\, where 1 0,<f>) = cos(0/2) |T)+ 
e 8 ^ sin(0/2) \l) indicates an inital spin-1/2 coherent state 
in the qubit basis (the two-dimensional Hilbert space 
spanned by the (1,1) charge state and spin-up on dot 1), 
U = H or T is the ideal intended single-qubit gate op- 

-i Jq* 3 dtri{t) is the true time 

evolution of the system under the time-dependent Hamil- 
tonian H(t). The overbar indicates a Gaussian average 
over fluctuations in the classical Zeeman field which 
reproduces the effects of hyperfine-induced decoherence 
due to an unknown static nuclear field when a ^ b2| 28 : 



eration, and U — Texp 



/(b 2 ) = J 



d 3 b 2 



(2tto~ ) 



3/2 



exp 



(b 2 - b 2 
2f7 2 



/(b 2 ). (6) 
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For a gated lateral GaAs quantum dot, a 

(b 2 — 2 = 3<7 2 due to hyperfme fluctuations has been 
measured, giving ctjv = 0.03 fieVM 
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Based on the above protocol for gating operations, and 
assuming a coherence time r c for the qubit spins, a suit- 
able parameter regime for high-fidelity single-qubit oper- 
ations is given by the following hierarchy: 



1/t c < 6f < 6| < t 12 < b\ < |e B | < |e A | 



(7) 



The first inequality in Eq. guarantees that a;- 
rotations are achieved with high fidelity at the operating 
point e = es- The second inequality allows for high- 
fidelity z-rotations at e = £a- The third and fourth in- 
equalities are required to ensure that b\ can be cancelled 
by exchange J ~ 2tf 2 /e, and the last two inequalities 
guarantee that the population of (0,2) (the double occu- 
pancy D ~ (ti2/e) ) remains small, which limits errors 
due to leakage and orbital dephasing (see below). When 
t c is dominated by hyperfine fluctuations, l/r c ~ <jn- In 
this case, we give a set of values for these parameters 
satisfying Eq. ([7]) in the caption of Fig. [3] The effec- 
tive Zeeman-field gradient given here could be achieved 
under the following circumstances: (a) a GaAs double 
quantum dot with the nuclei in dot 1 at near full po- 
larization, which would produce a maximum effective 
Zeeman splitting of b\ ~ 135 fj,eV (high polarizations 
could be achieved, e.g., through optical pumping 2 ^ or 
transport^), or (b) a nanomagnet neighboring a car- 
bon nanotube or graphene double quantum dot with g- 
factor g = 2 and inter-dot separation AL ~ 1 /im or an 
InAs nanowire double quantum dot with g-factor g = 8 
and inter-dot separation AL ~ 100 nm, either of which 
would require a magnetic-field gradient on the order of 
AB/AL ~ lT//im. Comparable field gradients have al- 
ready been achieved experimentally. 31 Alternatively, the 
ancillary spins could be polarized with the exchange field 
from a neighboring ferromagnet, high g-factor material, 
or stripline currents (see Fig. [I}. The values we have used 
for the detuning parameter e and tunnel coupling <i2 are 
of the same order as those used in previous experiments. 2 
Within the validity of the two-dimensional effective 
Hamiltonian U « exp {— iH c g (e)i s }, it is straightfor- 
ward but tedious to calculate rotation errors at e = £a 
(z-rotations: U = Rg(<f>)) and e = £b (^-rotations: 
U = Rx{4>)) using the expressions in Eqs. (|5I6[) . 3 - The 
error rate for z-rotations is dominated by the misalign- 
ment of the average field b 2 with the z-axis and is thus 
small in the ratio For a rotation by angle —<j> (to 

leading order in b'^/b^), this error rate is 



2 (b. 



sin — 



When A z =0 (at e = cb), the error rate for x-rotations 
is dominated by hyperfine fluctuations, and is therefore 
small in CT/v/fef. We find this error rate for an a;-rotation 
by angle —<j) (to leading order in V&f) i s 



(9) 
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Figure 3: (Color online) Error rates for H and T gates. 
For these plots we have chosen the parameters t±2 = 
100 ixeV, bf = 135 fieV, 6f = 10 fieV, 6| = 1/ieV, and 
eA ~ — 6meV. For the Hadamard gate, these values result in 
es = — £12/62 = ~ WOO fieV, a pulse time t p — 7rft/26| = Ins, 
and a rise/fall time r ~ 7rft/46| = 50 ps. Symbols give 
the results of numerical integration of the time-dependent 
Schrodinger equation for the Hadamard gate (H , circles) and 
7r/8 gate (T, squares). Lines give the estimates for gating 
error from Eq. (|10[) . 



We estimate the error in T using ry z ((/)) with <j> = ir/4. To 
estimate the error in H, we use Eq. ^ in combination 
with Eqs. @ and ([9]), assuming the errors incurred by 
each rotation are independent. These estimates give 



?f (f ) , (U = T), 

rf (f) +2^(1), (U = H). 



(10) 



From Eq. (jTUJ) we find the error rate for H is always 
larger than that for T and reaches a minimum at an 
optimal value of 6|. The optimal value of b\ and 77 at 
this point are: 



6 x,opt 



CM\a N , r,(br pt ) 



3 H\ 



(ii) 



(8) m 



where C is a numerical prefactor C — Wl/S + 7r 2 /48 ~ 
0.73. Using the measured value ctat = 0.03 fieV and 6| = 
10 [ieV, we find an optimized error rate of r\ ~ 10 -3 . Here 
we have included the most dominant error mechanisms. 
There are many other potential sources of error, which 
we discuss in the following. All numerical estimates are 
based on the parameter values given in the caption of Fig 



The error due to leakage to the (0,2) singlet state 
or misalignment of bi due to the hyperfine interac- 
tion in leading-order perturbation theory is given by 



(a N /bl) 2 , (t 12 /e A y 



10- 



If switching is done too slowly during the Hadamard 
gate, the qubit states will follow the adiabatic eigenbasis, 
introducing an additional source of error. We estimate 
this error to be 1 — P « a, where P = e~ a is the Landau- 
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Zeener tunneling probability, determined bys^ 



tt|6: 



i\ 2 



\dJ{t)/dt\ 



2tj 2 \Ae\ 



lo- 



in) 



Here, we have used dJ(t)/dt — 2ei 2 2 /e 2 3 , with |e| ~ 
|Ae|/r, where Ae = e^ — e B - In the opposite limit, 
a 3> 1, the qubit spin could be read out via charge 
measurements 2 - by sweeping slowly to large positive e, 
where the qubit state |t) would be adiabatically con- 
verted to the (0,2) ground-state singlet, or initialized by 
sweeping in the opposite direction (see Fig. Etc)). 

In systems with finite spin-orbit coupling, the 
transverse-spin decay time T 2 is limited by the energy 
relaxation time Ti (i.e., T 2 — 2T^ 3 ), so it is sufficient 
to analyze this error in terms of T\ . T\ in quantum dots 
can now be measured, 34 giving T\ ~ 1 ms at fields of 
B « 6T (g*^ B B ~ 135 fieV)M This value gives an error 
estimate on the order of t s /Ti < 10 -6 for a switching 
time i s — 1 ns. 

Finally, rapid voltage-controlled gating in this scheme 
is made possible only because the electron spin states 
are associated with different orbital wave functions dur- 
ing pulsing, which also makes these states susceptible 
to orbital dephasing. The associated dephasing time 
is, however, strongly suppressed in the limit where the 
double occupancy is small: D as (ti2/e) 2 <C 1. In par- 
ticular, the dephasing time for the two-electron system 

is > D -2 ^ 1 ^,— where at Ins 36 is the single- 
electron dephasing time in a double quantum dot. This 
gives an error estimate of t s /r^ < 10~ 4 , using t s as 1 ns 
and D ~ 10~ 2 at the operating point e = e B . It should be 
possible to further suppress orbital dephasing by choos- 
ing the operating point e B to coincide with a "sweet 
spot", where dJ(e B )/de = O^^l 

Numerical analysis-To confirm the validity of the ap- 
proximations made here and to verify the smallness 
of error mechanisms associated with leakage and finite 
pulse times, we have numerically integrated the time- 
dependent Schrodinger equation for the Hamiltonian 



given in Eq. ([T]) in the basis of the (0,2) singlet state 
and four (1,1) states (including spin). We have used the 
pulse scheme described following Eq. and evaluated 
the gate error rates for T and H from the fidelity in Eq. 
(O . For the Hadamard gate we used the symmetric pulse 
shape 



e(t) 



eo 



4^ tanh ( 



2\t-2r] 



,0 < t < I 



eo 



^tanh( 2 ' t -- 2T - t n ,^<t<t ! 



(13) 



where eo = (e^ + eg) /2 and Ae = e B — £A- The pulse 
time t p and rise/fall time r = (t s — t p )/4 were fixed using 



7T 



Mt)dt, 



(14) 



where the solution to the above integral equation was 
found numerically. The results of our numerics are shown 
in Fig. [3] To implement the integral (Eq. ^) numeri- 
cally, we have performed a Monte Carlo average over 100 
Overhauser fields, sampled from a uniform Gaussian dis- 
tribution using the experimental value er/v = 0.03 /ieV. 
Error bars due to the finite sample of Overhauser fields 
are smaller than the symbol size. We find good agree- 
ment between the analytical and predicted error rate for 
T in the limit of large 6f (the saturation value for rj at low 
6f is consistent with our estimates of ~ 10~ 4 for error due 
to leakage). Additionally, we find reasonable agreement 
with our estimate for the -ff-gate error rate, confirming 
that we have identified the dominant error mechanisms. 
This gives us confidence that an error rate on the order of 
~ 10~ 3 should be achievable with this proposed scheme. 

Acknowledgments-We thank G. Burkard, H.-A. Engel, 
M. Friesen, V. N. Golovach, J. Lehmann, D. Lidar, B. 
Trauzettel, and L. M. K. Vandersypen for useful dis- 
cussions. We acknowledge financial support from JST 
ICORP, EU NoE MAGMANet, the NCCR Nanoscience, 
and the Swiss NSF. 



1 D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 
(1998). 

2 J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, 
A. Yacoby, M. D. Lukin, C. M. Marcus, M. P. Hanson, 
and A. C. Gossard, Science 309, 2180 (2005). 

3 F. H. L. Koppens, C. Buizert, K. J. Tielrooij, I. T. Vink, 
K. C. Nowack, T. Meunier, L. P. Kouwenhoven, and 
L. M. K. Vandersypen, Nature 442, 766 (2006). 

4 A. M. Steane, Phys. Rev. A 68, 042322 (2003). 

5 E. Knill, Nature 434, 39 (2005). 

6 H.-A. Engel and D. Loss, Phys. Rev. Lett. 86, 4648 (2001). 

7 R. Hanson and G. Burkard, Phys. Rev. Lett. 98, 050502 
(2007). 

8 J. Kyriakidis and G. Burkard, Phys. Rev. B 75, 115324 
(2007). 



9 E. I. Rashba and A. L. Efros, Phys. Rev. Lett. 91, 126405 
(2003). 

D. Stepanenko and N. E. Bonesteel, Phys. Rev. Lett. 93, 
140501 (2004). 

1 C. Flindt, A. S. S0rensen, and K. Flensberg, Phys. Rev. 
Lett. 97, 240501 (2006). 

2 V. N. Golovach, M. Borhani, and D. Loss, Phys. Rev. B 
74, 165319 (2006). 

3 Y. Tokura, W. G. van der Wiel, T. Obata, and S. Tarucha, 
Phys. Rev. Lett. 96, 047202 (2006). 

4 L.-A. Wu, D. A. Lidar, and M. Friesen, Phys. Rev. Lett. 
93, 030501 (2004). 

5 W. A. Coish and D. Loss, Phys. Rev. B 72, 125337 (2005). 

6 X. Hu and S. Das Sarma, Phys. Rev. Lett. 96, 100501 
(2006). 



■5 



17 M. Stopa and C. M. Marcus, arXiv.orgxond-mat/0604008 
(2006). 

18 M. Friesen, P. Rugheimer, D. Savage, M. Lagally, 
D. van der Weide, R. Joynt, and M. Eriksson, Phys. Rev. 
B 67, 121301 (2003). 

19 N. Mason, M. J. Biercuk, and C. M. Marcus, Science 303, 
655 (2004). 

20 M. R. Graber, W. A. Coish, C. Hoffmann, M. Weiss, 
J. Furer, S. Oberholzer, D. Loss, and C. Schonenberger, 
Phys. Rev. B 74, 075427 (2006). 

21 B. Trauzettel, D. V. Bulaev, D. Loss, and G. Burkard, Nat. 
Phys. 3, 192 (2007). 

22 F. H. L. Koppens, J. A. Folk, J. M. Elzerman, R. Han- 
son, L. H. W. van Beveren, I. T. Vink, H. P. Tranitz, 
W. Wegscheider, L. P. Kouwenhoven, and L. M. K. Van- 
dersypen, Science 309, 1346 (2005). 

23 W. G. van der Wiel, S. de Franceschi, J. M. Elzerman, 
T. Fujisawa, S. Tarucha, and L. P. Kouwenhoven, Rev. 
Mod. Phys. 75, 1 (2003). 

24 W. A. Coish and D. Loss, arXiv:cond-mat/0606550 (2006), 
to appear in the Handbook of Magnetism and Ad- 
vanced Magnetic Materials, vol. 5, eds. S. Parkin and D. 
Awschalom, Wiley. 

25 G. Burkard, D. Loss, and D. P. DiVincenzo, Phys. Rev. B 
59, 2070 (1999). 

26 E. A. Laird, J. R. Petta, A. C. Johnson, C. M. Marcus, 
A. Yacoby, M. P. Hanson, and A. C. Gossard, Phys. Rev. 



Lett. 97, 056801 (2006). 

M. A. Nielsen and I. L. Chuang, Quantum Computation 

and Quantum Information (Cambridge, 2001). 

W. A. Coish and D. Loss, Phys. Rev. B 70, 195340 (2004). 

A. Imamoglu, E. Knill, L. Tian, and P. Zoller, Phys. Rev. 
Lett. 91, 017402 (2003). 

M. S. Rudner and L. S. Levitov, arXiv.org:cond- 
mat/0609409 (2006). 

J. Wrobel, T. Dietl, A. Lusakowski, G. Grabecki, K. Fronc, 
R. Hey, K. H. Ploog, and H. Shtrikman, Phys. Rev. Lett. 
93, 246601 (2004). 

C. Zener, Proc. Roy. Soc. Lond. 137, 696 (1932). 

V. N. Golovach, A. Khaetskii, and D. Loss, Phys. Rev. 

Lett. 93, 016601 (2004). 

J. M. Elzerman, R. Hanson, L. H. Willems van Beveren, 

B. Witkamp, L. M. K. Vandersypen, and L. P. Kouwen- 
hoven, Nature (London) 430, 431 (2004). 

S. Amasha, K. MacLean, I. Radu, D. Zumbuhl, M. Kast- 
ner, M. Hanson, and A. Gossard, arXiv:cond-mat/0607110 
(2006). 

T. Hayashi, T. Fujisawa, H. D. Cheong, Y. H. Jeong, and 
Y. Hirayama, Phys. Rev. Lett. 91, 226804 (2003). 
The integral over initial states in Eq. (5) can be replaced by 
a discrite sum, which presents a significant simplification 
(see M. D. Bowdrey, D. K. L. Oi, A. J. Short, K. Banaszek, 
and J. A. Jones, Phys. Lett. A 294, 258 (2002)) 



